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We derive corrections to the JIMWLK equation in the regime where the charge density in the 
hadronic wave function is small. We show that the framework of the JIMWLK equation has to be 
significantly modified at small densities in order to properly account for the noncommutativity of the 
charge density operators. In particular the weight function for the calculation of averages can not be 
. real, but is shown to contain the Wess-Zumino term. The corrections to the kernel of the JIMWLK 

evolution which are leading at small density are resummed into a path ordered exponential of the 
functional derivative with respect to the charge density operator, thus hinting at intriguing duality 
£S| . between the high and the low density regimes. 
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I. INTRODUCTION. 



The recent years have seen vigorous activity in an attempt to understand high energy scattering in QCD from first 
principles. The precursors of this recent wave are the classic papers Q and Q as well as later work and 0. 
■ The equation that governs the evolution of the S-matrix of a small projectile with energy was derived by BalitskyQ, 
and its mean field version by Kovchegov|y|. This equation takes into account perturbative evolution of the projectile 
wave function supplemented by multiple scatterings of the projectile partons on the target. A lot of numerical and 
ON ■ analytic work in recent years has been carried out to understand the properties of the solutions of this equation both 
in the ultraviolet and infrared. The dependence of the saturation momentum on rapidity 0, the property of geometric 
scaling 9j , the power growth of the total cross section § and the disappearance of the Cronin effect with energy [To| 
1^ | are some stark examples of the physical information obtained from this equation so far. 

A complementary approach to the high energy scattering problem was initiated in Il2| . It yields the nonlinear 
evolution equation (JIMWLK equation) that governs the change of the correlation functions of the color charge density 
in the hadron wave function. The parameter of the evolution, as before is the rapidity of the hadron. In this approach 
II 1 it is assumed that the charge density in the wave function is high. It has been shown in the last reference in [TJ and 
Q_i in ^K see a l so E3) that when understood as the equation for the wave function of the large target, for processes 
when this target is probed by a small perturbative projectile, the JIMWLK equation is equivalent to that of Balitsky 
as far as the calculation of S'-matrix is concerned. 

One of the original motivations to consider the high energy scattering, is the question how does the scattering 
amplitude approach the unitarity limit. The Balitsky-Kovchegov and JIMWLK equations indeed lead to unitary 
amplitude: at large energies the S'-matrix approaches zero (not to be confused with the unitarity of the total cross 
section in the sense of the Froissart bound, which is not achieved in this framework Q). However it is clear that the 
way the amplitude approaches the unitarity limit is not properly described by these equations. In the framework of 
the Balitsky-Kovchegov the reason is that the evolution of the projectile wave function is linear, that is gluons are 
emitted independently by the "valence" partons in the projectile. The unitarization of the amplitude is then achieved 
only due to the multiple scatterings on the target. However it is clear that when evolved to large enough rapidity, the 
wave function of the projectile becomes dense, and at that point the gluons have to be emitted coherently from the 
charge density distribution in the target, rather than independently from every parton. At high density the gluons are 
produced less efficiently due to these wave function saturation effects, and this furnishes an additional mechanism for 
the unitarization of the scattering amplitude over and above multiple scatterings. When the target is large and the 
projectile is small, the multiple scattering mechanism is more efficient. However in this situation one already at low 
energy starts from the almost vanishing scattering matrix (black target), and so can not sensibly study the approach 
to the black limit. On the other hand, for a small target, the unitarity corrections due to the wave function saturation 
effects in the projectile become important at the same rapidity as those due to multiple scatterings, and the approach 
to the unitarity limit is described incorrectly when they are omitted. 

In the framework of the JIMWLK equation one is faced with a similar problem. The equation is derived in the 
limit of high density. Although its low density limit reduces to the linear BFKL equation, there is no reason to believe 
that it describes correctly the evolution in the intermediate regime, where the charge density is neither small nor 
large. Therefore one cannot consistently use the JIMWLK equation to follow the evolution of the wave function of a 
" dilute" system all the way until it becomes " dense" . 

The main physical effects that take place in the intermediate density regime are due to the so-called Pomeron loops. 
Those are the processes where the gluons, which are emitted in the wave function at an earlier stage of the evolution, 
subsequently " drop out" of the evolution as their color is " bleached" by other gluons and they cease their ephemeral 
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existense. The JIMWLK equation contains only part of the Pomeron loops. Although both, the gluon emission and 
the gluon disappearance processes have their place in the JIMWLK evolution, they are not described properly at small 
and intermediate densities (rapidities) . Proper inclusion of these processes must also restore the t-channel unitarity 
which is not preserved by the JIMWLk equation |Lj|. One would like to have a better handle on the Pomeron loops, 
since they are clearly important for the evolution of the scattering amplitude. 

To that end one has to understand the evolution of the hadronic wave function at arbitrary color charge den- 
sity. This is a very challenging problem which at present is not solved. Recently a first step in this direction was 
undertakenfill . [lq . [l7j . These works derive a correction to the JIMWLK equation, which can be described as the 
first term in expansion in powers of the functional derivative with respect to the color charge density. In the present 
paper we go beyond this first step and resum an infinite number of such terms. The terms we find do not give the 
complete kernel of the JIMWLK equation at arbitrary density. We do however resum all terms which are leading in 
the low density regime (plus a little bit extra - see below). Along the way we find that the framework of the JIMWLK 
equation has to be modified to accomodate the low density regime. In particular, since one cannot treat the charge 
density operators as commuting when the density is not large, the quantum averages can not be written simply as 
weighted averages, (functional integrals) over the real functional JF"[p(a;)]. Rather the charge density variable in the 
functional integral must be considered as dependent on an extra variable t and the weight functional W must contain 
a phase which is given by the Wess-Zumino term. We also show that this auxiliary variable t plays exactly the same 
role as the longitudinal coordinate x~ in |18| and |l2j . namely it resolves the ambiguity associated with the ordering 
of quantum fields. Thus for all intents and purposes t can be identified with x~ . 

The correction terms we derive are resummed in a compact and suggestive expression. It has a very interesting 
structure which hints at an intriguing duality between the low and high density regime. 

This paper is organized as follows. In section 2 we derive the correction to the correlation functions of the charge 
density operator in terms of averages of quantum operators. In section 3 we translate this into the corrections to the 
functional evolution equation for the functional W[/a]. We also explain how the extra variable t and the Wess-Zumino 
term arise in this context. Finally in section 4 we discuss several issues and questions which are prompted by our 
results. 



II. HIGH ENERGY HADRONIC WAVE FUNCTION AND THE EVOLUTION OF THE 
CORRELATORS OF COLOR CHARGE DENSITY. 

The JIMWLK evolution equation is a functional evolution equation for the weight functional W[p a (x)}. It was 
shown in that the form of the evolution is somewhat simpler when W is considered to be a function of the phase 
of the scattering matrix a rather than p. For our purposes in this paper however we find it more convenient and 
straightforward to use the original formulation in terms of The functional W has the meaning of the probability 

density to find a given configuration of the charge density, so that expectation value of any observable 0[p] in the 
hadronic wave function is given by 

(O) = J Dp a O[p{x)]W[p{x% J Dp a (x)W[p(x)} = 1 (2.1) 

Here Xi, i = 1,2 are the transverse coordinates, and p a {x) is the surface charge density, defined as the integral of the 
three dimensional charge density over the longitudinal extent of the hadron 

p a {x) = J dx-p a (x,x-) (2.2) 

The JIMWLK evolution equation reads 

d :W[p] =aj [ # x #yx*\p]**— + [ d 2 xa a (x)-^\ W[p] (2.3) 



dY irs s \J y/V irs p a {x) 5p b (y) J y 'p a (x) 

The explit form of the functionals x an( i a 1S n °t important for our present purposes, except for the fact that they 
both have expansion in powers of p, so that at small p one has x = 0(p 2 ) an( i a = O(p). 

The equation was derived for parametrically large densities p oc — . We will now extend the derivation to include 
also the small density region p oc 1. The easiest way to do this is to consider directly the evolution of the hadronic wave 
function. Although the original derivation in |Tl| was given in the language of the path integral, it can be formulated 
directly in terms of the wave function [l9| . One starts at the initial rapidity with a wave function \^{Y)) — \v) which 
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contains only "valence" degrees of freedom, namely those with the longitudinal momentum k + above some cutoff 
value A. This wave function defines the charge density correlation functions 



(v\p ai (x 1 )...p a "(x n )\v) 



(2.4) 



When boosted by a small amount, the valence wave function gets dressed by a cloud of the Weiszacker- Williams 
gluons. The change of the wave function can be calculated explicitly. The evolved wave function has the following 
structure 113 



\V(Y + 5Y)) 



/,/-,(///(.-)//;(.-)) 



i / d 2 xbf(x) 



dk+ 



(1- S Y)A 7rV2|frf|l/2 



al a (k+,x)\ B(ala)\v) (2.5) 



Here the creation operators a^(fc + ) create gluons with soft momenta, which are not present in the valence state \v) >. 
The field b depends only on the valence degrees of freedom. It is determined as the solution of the " classical" equation 
of motion 



dibt + ge^ixMx) = gp a v (x) 



Am 



t ahc b h rt 







(2.6) 



This is precisely the "classical background field" that appears in 11]. The only subtlety is that the commutator 
in the first equation was dropped in ^l]], which is indeed appropriate for large charge densities. In this case it is 
straightforward to see that the commutator term is 0(a s ) correction to the first term in the equation and can therefore 
be nerglected. In fact for small charge densities, p — O(l) the field b is of order g, and so the commutator term is 
again negligible. We keep it here for completeness, as it can in principle be important in the intermediate regime. 
The p v in the right hand side of eq. I|2.(j[) is the color charge density of the valence degrees of freedom only, that is of 
the gluons with longitudinal momentum above the cutoff A. We stress again, that in the present paper we consider 
p v and bi to be fully quantum operators which act on the Hilbert space of the valence degrees of freedom. 

The term in the curly brackets in ca. H2.5JI is simply an expansion to first order in SY of the coherent operator 



C = exp • 



d 2 xb a l (x) 



dk+ 



al a (k+,x)+aU^ 



{ 1-8Y)A 7^+11/2 



which affects the shift of the soft modes of the gluon field by the "classical field" bf(x) 

&A$(k+,x)C = A1{k+,x) + pr&ftx) 



(2.7) 



(2.8) 



In ea. H2.5[) we have only kept C to first order in SY, since we mean to calculate the first derivative with respect to 
Y, and thus the higher order terms do not contribute. The last ingredient in ea. (|2.5(l that requires explanation is the 
operator B. This operator is the exponential of the quadratic form of a'a and a''s which also depends on the valence 
charge density p v . It's role is to perform the Bogolyubov transformation on the soft gluon fields, mixing the creation 
and annihilation operators |_19| . The two operators in ea. l|2.5[) . C and B correspond precisely to the two steps in the 
calculation of pd| - the expansion of the quantum fields around the classical background b, and subsequent integration 
over the soft modes keeping only quadratic terms in the Lagrangian. The more complicated mathematics of |Tl| (the 
one loop integration over the soft modes) is coded in the operator B. However luckily for us, it has the property 
that for small densities it becomes a unit operator B(p = 0) = 1. Thus in the small density regime it only brings in 
perturbative corrections, and therefore in the following we set it to unity. 

In fact, once we set B = 1 the wave function ea. (|2.5|l becomes precisely the one used to derive the Balitsky equation 
in pcf. with the substitution of the lowest order Weiszacker- Williams field J d 2 y ^JZyyi P a (y) by bf(x). Of course, b 
which solves the nonlinear equations ea. (|2.6|l . when expanded to first order in p, reduces to this expression. 

Given the evolution of the wave function, we can now calculate the evolution of the correlation functions of the 
charge density. Since the boost operation "opens up" the Hilbert space of the soft modes, we have to consider also 
the contribution of these modes to the charge density operator 



P a (x)=p a v (x) 



dk+af(k+,x)T b a c a c l (k+,x) 



(1-5Y)A 



(2.9) 



where T h a c = if abc is the SU(N) generator in the adjoint representation. 
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Sandwiching the n-th power of p in the wave function eq. 1)2. 5f) and differentiating with respect to Y we obtain 

±(p^( Xl ),..p^( Xn )) = i{l [ d 2 x([b?( X ),{bUx),P ai (xi)-P a "(Xn)})+ (2.10) 



J d 2 xY,8 2 {x ~ Xi)S 2 (x - Xj ){T a <T a i) bc 



x(&^)^i)---P a! -^-i)p a ^^^ 



+ J d 2 x6 2 (x-x l )...S 2 (x-x n )(T a \..T a -) bc (b b (x)b c t (x))} 

where now the averages are over the valence state \v) and the charge density operators are those of the valence degrees 
of freedom only. 

The i-th line in ea. (|2.11fl contains terms with i operators p deleted and substituted by the single gluon charge 
density operators S 2 (x — Xi)T ai . This structure is completely transparent. One step in the evolution corresponds to 
emission of one gluon. The first term in ea. (|2.11|l corresponds to the rotation of the valence charge density by the 
charge of the emitted gluon. All other terms are due to the direct contribution of the color charge density of the 
emitted gluon itself. 

Note that the first (double commutator) term on the right hand side of ea. l|2.11|) despite appearances starts with 
the same power of p as the left hand side. Both, in the weak and strong field limit the relevant commutator can be 
written as 

lb-(x),p b (y)] = j^f cbd p d (y) (2.ii) 



and 

<«>i(s) 
6p c (y) 



= f) 



Di 



d ■ D 



(x,y) (2.12) 



Thus each commutator reduces the lowest power of p by one. This term is contained in the JIMWLK evolution at 
large fields, as it can be written as the contribution to the kernel with two functional derivatives with respect to p 
multiplied by a function of p whose Taylor expansion starts at p 2 . It will therefore not interest us in the following, 
and we will concentrate our attention on the rest of the terms in ea. l|2.11|l . 

These terms have insertions of T a in place of p a , and so one naturally would like to write them down as functional 
derivatives with respect to p. Here however one faces the problem, that the charge density operators do not commute 
with each other, and thus the ordering of p's and 6's in ea. (|2.11|l is important. Fortunately there is a known way of 
representing correlation functions of noncommuting variables as a functional integral. We describe it, as well as its 
application to the present case in the next section. 



III. THE FUNCTIONAL INTEGRAL REPRESENTATION AND THE WESS-ZUMINO TERM. 



The general problem we have to address is how to represent the correlators of the SU(N) generators p a (x) in terms 
of the functional integral. We will describe this construction in the case of SU(2). Let us first forget about the x 
dependence and also assume that the gen erators is in the fixed representation of spin J. The construction for this 
case has been worked out in detail in [21j . It is is based on the observation that instead of considering the ordered 
product of the generators p a in the representation J, one can consider the correlation function 



(p a (h)p b (t 2 ) . . . p c (t k )) — J k (n a (h)n b (t 2 ) . . . n c (t k )) 



(3.1) 



Dn{t)n a (t 1 )n b (t 2 ) ■ ■ - n c (t k ) exp 



4 I d A £,e aP e abc n a d a n b dp 



where J is the spin of representaion (i.e. for fundamental representation J = 1/2), n a (t) is a unit vector n a n a = 1 
living on a contour C (t is a coordinate on the contour) and S is an arbitrary two-dimensional surface with the 
boundary C = <5S. Despite the appearance, the two-dimensional Wess-Zumino action 



S[n] 



d 2 Ze a0 e abc n a d a n b dpn c 



(3.2) 
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depends only on values n a (t) at the boundary and not on the values n takes on the surface E. 

The coordinates U on the left hand side of eg. (|3.2I) are only important as indicators of the ordering of the operators. 
Thus the correlation function of n's must only depend on the ordering of the time coordinates, and not on their values. 
This is indeed the case. To see this note that the variation of the Wess-Zumino action under arbitrary transformation 
of the fields n is 



<)S = <j> dte abc n a d t n b 5n~ 



Now let us perform an infinitesimal SU{2) transformation 

n a (t) -> n a {t) + e abc fl b (t)n c (t) 



(3.3) 



(3.4) 



The action change under this transformation is SS = — § c dth a (t)VL a (t) Performing this transformation in the path 
integral for the correlation function eq. l|3.2|l we obtain 



J^(n a (t)n 6 (tx) . . . n c {t k )) = - (n f (t)n b (h) . ..n d (U).. . n c (t k )) 



(3.5) 



i=l 



where n d (ti) means the exclusion of this term from the products of the fields in a correlator. This establishes the 
piecewise constant nature of the correlation function. Moreover, remembering that in the equal time limit the time 
derivative of the correlation function reduces to the equal time commutator (as the path integral represents T-ordered 
products), one concludes immediately from ea. l|3.5|l that the following equal time commutaion relations hold 



J \n ,n = le n 



This establishes that the following identification is valid 



P 



Jn a 



(3.6) 



(3.7) 



The fact that the charge density depends on the transverse coordinate simply means that we have to make the 
unit field n a also x-dependent, but the Wess-Zumino term is strictly local in x (since otherwise we would introduce 
noncommutativity between p's at different points in the transverse plain). Finally we should also allow to consider 
states with different representation of the SU (2) group. This is achieved by allowing J to be distributed with some 
arbitrary weight and supplementing the integration measure over n a by summation over all half integer J. 

All said and done we see that we can represent the calculation of charge density correlators by the following 
functional integral 



{p a ^{x 1 )...p a -{x n )) = / dp{x,t)p a ^x 1 ,t 1 )...p a "(x n ,t n )W[p] 



with the measure of integration being understood as 



Jdp(x,t)... = J ^j( x )U Xt tdn(x,t)ex.p i J d 2 x^^- J d 2 ^e a 



^e abc n a {x, t)d a n b (x, t)d p n c (x, t) 



(3.8) 



(3.9) 



The generalization of this procedure to the SU (N) group is somewhat cumbersome and we will not describe it 
in detail, but rather only indicate the flow of the argument. It goes along the following lines. The Wess-Zumino 
term in eq. (|3.2|) can be thought of as the flux of the magnetic field of the 'tHooft-Polyakov monopole through the 
contour E. This monopole would sit in some three dimenisonal space, two of whose dimensions are spanned by the 
coordinates on the surface E. It is of course purely fictitious and should not be given any physical significance beoynd 
a simple mnemonic to describe the mathematical structure of the relevant Wess-Zumino term. The conservation of the 
magnetic flux of this monopole modulo 47r is the reason why the Wess-Zumina term does not depend on the surface 
E, but rather only on its boundary, as long as the coefficient J takes half integer values. In the SU(N) group there 
are N — 1 such independent magnetic monopoles, whose magnetic fields and therefore fluxes can be expressed in terms 
of the adjoint "Higgs field" of unit length n a , a — 1, n 2 — l j22|. Thus there are N — 1 independent Wess-Zumino 
terms that can be written for the SU(N) group. The coefficients of all these terms have to be quantized for the same 
reason as J, and this quantization corresponds to the quantization of the eigenvalues of the generators of the Cartan 
subalgebra of the SU (N) group. To handle the generators of SU (N) in a given representation one again has to endow 
p a with a dependence on one extra coordinate t, and introduce N — 1 Wess-Zumino terms in the integration measure. 
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The coefficients of these terms should be equal to the values of the Cartan generators in the highest weight state of 
the given representation. Taking into account x dependence and allowing for the variation of the representation is 
achieved in the same way as in ea. (|3.9|) . 

We note that in the limit of large representations, J — > co, the Wess-Zumino term imposes the constraint h a {t) = 0. 
Thus the ^-dependence of the unit vector n is frozen, n a (t) = n a , so that one recovers the formulation of eq. l|2.1[l . 

Practically speaking therefore we see that taking into account noncommutativity of p leads to two major changes, 
first the field in the functional integral becomes dependent on extra coordinate t, and second the measure of the 
integration becomes complex with the phase given by the Wess - Zumino term. 

Nevertheless, given the functional integral representation, we can conveniently rewrite the evolution equations 
ea. (|2.11l) as a functional derivative operator acting on W [/£>]. It is in fact easy to see that the whole hierarchy of 
cos. <|2.11|) can be written in a very coincise form 



d 2 xb b (x 1 t^ 



^)T£(b»(x)p a i(x 1 )...p a ^(x^ 1 )p a ^(x t+1 )...p a "(x n mx)) = 



(3.10) 



-oo)T b a c ^(:M ^ +oc) / dt 



5p a (x,t) 



{p ai (x 1 ,t 1 )...p a "(x n ,t n )}); 



f d 2 x^2S 2 (x-x l )S 2 (x- Xj )(T a *T a >) bc 
J i<j 



x{b\{x)p a -{xx)... P a ^{ Xi -x)p a ^{x i+ i)...p^ 



( / sfcjfot -cx))(T a T c V& l c (a;,< -» +oo) 
etc... 



dt 



-foe 



dt' 



8p a (x,t) Sp b (x,t') 



-{p ai (x u t 1 )...p a " {x n ,t n )}) 



Thus somewhat surprisingly we can rewrite the right hand side of eg . (|2. 1 lfl (subtracting the first term) as 



( / d 2 xb b (x,t^> -co) 



Pexp 



dtT a - 



Sp a (x,t) 



bl{x,t -» +oo){ P a ^xi,t 1 )...p a -(x n ,t n )}) (3.11) 



where P denotes the path ordering along t. Note that the linear and quadratic terms in the expansion of the path 
ordered exponential are contained already in the JIMWLK equation (modulo the noncommutativity of p's). The 
rest of the terms are subleading at large p. However at p = 0(1) all the terms in ea. (|3.11|l are of the same order 
and should be kept. Also note, that by keeping bf as the full solution of the equation ea l2.6[) and not its leading 
perturbative term, we are resumming some terms that are subleading in the low density limit. Partially integrating 
the functional derivatives in the functional integral representation ea. l|3.8[l we can finally rewrite the evolution as the 
functional equation for W 



d_ 

w 

1 

H — 



W[p(x,t)] = a s \ Ojimwlk( P , ^)W[p(x,t)] 



d 2 xb\(x, t — » — co) 



P exp < — 



dtT a 



5p a (x,t) 



- 1 



b\(x,t 



W[p(x,t)} (3.12) 



where P is the path ordered exponential with the linear and quadratic terms subtracted. 
This is the main result of the present paper. 



IV. DISCUSSION. 



There are several interesting questions that arise from the previous derivation. First, does the variable t has a 
physical meaning beoynd being a useful tool to represent correlators of noncommuting variables. We believe that the 
answer to this question is in the affirmative. We remind the reader that one already had a need in introducing path 
ordering in the discussion of high energy evolution. In particular it was realized in 18] that this was necessary in order 
to solve the classical equations of motion for the gluon field in the presence of the charge density p. The equation 
considered in ^3 was identical to ea. H2.(j|) . except b was considered to be a classical field which depended on x as 
well as on the longitudinal coordinate x~ . It then turned out that only one contribution to the commutator term in 
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the first equation of ea. l|2.6f) was important. This contribution was the one where one of the 6's had the longitudinal 
coordinate slightly smaller than that of the other b. The equation then was formally solved as 

bi(x, x~) = iW(x, x~)diU(x, x~), U{x, x~) — Pexp{i / dx~a(x, x~)} (4-1) 

J — oo 

with the function a. itself defined in terms of the matrix U and the charge density p as a = [UpU^]. The important 
element in these expressions is the path ordering with respect to x~ . Now returning to our present framework we can 
ask ourselves what is the solution of the operator equation eq. i|2.6[l that appears in the functional integral expression 
ea. (|3.12|l . It is obvious, that in order to solve this equation one has to follow the procedure identical to that in [Isj . 
One has to endow b and p with the extra coordinate t, the commutator term then has the t coordinate of the first 
6-factor slightly smaller than that of the second factor, as the bordering in the classical equation simply reflects the 
operator ordering in the operator equation. From that point on the solution is identical to that of |18| with the path 
ordering in x~ replaced by the path ordering in t. 

We conclude therefore that the ordering variable t that we have introduced in this paper is identical to the longi- 
tudinal coordinate x~[23|. 

One than immediately is tempted to ask: what about the other variable which is introduced to define the Wess- 
Zumino term? If t is x~ , then the other one may be x + . The connection here is more tenuous, but we believe 
it is true. The Wess-Zumino term can be thought of as the Berry phase [24j for the state in the J-representation 
of the charge density operator. The state depends on x~ as a parameter, and its phase changes as this parameter 
changes adiabatically along the hadron. The Berry phase arises as a topological part of the time integral of the 
simplectic form in the action. The QCD functional integral in the light cone gauge (to which the functional integral 
ea. (|3.2|) is supposed to be an approximation) is in fact a phase space path integral, and contains the simplectic form 
F +l F~ l = d + bid~bi. It is likely that the Wess-Zumino term in ea. l|3.2|l should be properly understood as the time 
integral of this simplectic form, and thus the second coordinate which appears in its definition is indeed the time 
variable x + . This is an interesting question, which should be understood better. 

Another natural question that arises, is whether one can somehow avoid introducing the noncommuting variables 
alltogether, and thus get rid of the awkward phase in the "weight functional" W . Unfortunately we do not believe this 
possible. Suppose we are interested in the expectation values of some set of observables {Oi}. If all these observables 
are mutually commuting, one can choose a basis in the Hilbert space (spanned by a set of coordinates which commute 
with {Oi}) such that all the expectation values are given as averages over the real measure. The measure is the square 
of the wave function in the basis we have chosen. However if we are also interested in averages of other variables 
which do not commute with Oi, the calculation of the expectation value will necessarily probe also the phase of 
the wave function. The question is thus whether in high energy QCD we are only interested in mutually commuting 
observables. The observables we would like to calculate in the hadronic wave function are the averages of the S'-matrix 
of fast particles scattering on it. In general these observables do not commute. Physically it simply means that the 
probability for scattering of two particles depends on the order in which these particles scatter, since the first one to 
scatter perturbs the target fields. When the number of the projectile particles is small and the target fields are large, 
this small perturbation is a subleading effect and thus can be neglected. However for targets which are not dense this 
noncommutativity is important, and so the phase of the wave function comes into play. 

Perhaps the most interesting question that is suggested by our result eq. (|3.12[1 is whether there is some well defined 
duality transformation that maps the high density regime into the low density regime. Recall, that in the high density 
JIMWLK expression the prominent element is the matrix U defined as 

U = Pexp{i J dx~a{x,x^)} (4.2) 

The relation between a and p is fairly complicated, but at low density it simplifies considerably 

g 2 

a(x,x ) =io W density Tp(x,y)p(y,x ) (4.3) 

Thus it looks like the path ordered exponential of density and the path ordered exponential of functional derivatives 
with respect to density roughly exchange their roles in the high and low density regimes. Of course the swap is not 
so starightforward, but nevertheless the appearance of the path order exponential of j- in eq. H3.12fl is very suggestive 
that such a transformation can be defined. One could hope that the knowledge of this transformation will help find 
the exact expression for the operator of the high energy evolution valid at arbitrary density p. 

Finally we comment on the relation of our results to those presented in [Tt| . The paper is not explicit about 
the operator ordering of the operators involved, and that makes the detailed comparison a little ambiguous. Also, the 
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derivation in [16| are only given in the large N c approximation which relies on the dipole model. However it is easy 
to see that the general structure of their result is the same as our expression ea. (|3.12|) expanded to fourth order in 
the functional derivative. To establish this one has to use the perturbative relation between a and p ea. (|4.3[) and the 
commutator of the field b with p ea. (12. 11(1 taken to lowest order in the strong coupling constant. Additionally one 
has to assume that the funtional W depends only on the "dipole cross section" s(x,y) = -ktr[U'(x)U(y)]. 

Acknowledgments 

We thank Urs Wiedemann for many discussions on subjects directly relevant to this work. 



[1] L.V. Gribov, E. Levin and M. Ryskin, Phys. Rep 100:1,1983; A. Mueller and J.W. Qiu, Nucl. Phys. B286:427 (1986). 
[2] E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977), 199; Ya.Ya. Balitsky and L.N. Lipatov, Sov. J. 
Nucl. Phys. 28 (1978), 22. 

[3] A. Mueller, Nucl.Phys. B335:115 (1990); Nucl. Phys. B415:373-385,1994 ; Nucl. Phys. B437:107-126,1995 
[4] L. McLerran and R. Venugopalan, Phys. Rev. D49:2233-2241,1994 ; Phys.Rev.D49:3352-3355,1994 
[5] I. Balitsky, Nucl. Phys. B463:99,1996; Phys. Rev. Lett. 81:2024-2027,1998; Phys. Rev. D60:014020,1999. 
[6] Y. Kovchegov, Phys. Rev. D60:034008,1999. 

[7] E. Levin and K. Tuchin; Nucl.Phys. B573 (2000) 833-852; E. Levin, M. Lublinsky, Nucl.Phys. A696 (2001) 833-850 
[8] A. Kovner and U. Wiedemann, Phys.Rev. D66 (2002) 051502; Phys.Rev. D66 (2002) 034031; Phys.Lett. B551 (2003) 
311-316, 

[9] E. Iancu, K. Itakura and L. McLerran; Nucl.Phys. A708 (2002) 327-352; 
[10] J. Albacete et.al, Phys.Rev.Lett. 92 (2004) 082001; D. Kharzeev, Yu. Kovchegov and K. Tuchin; Phys.Rev. D68 (2003) 
094013 

[11] J. Jalilian Marian, A. Kovner, A.Leonidov and H. Weigert, Nucl. Phys.B504 415 (1997); Phys. Rev. D59 014014 (1999); 
J. Jalilian Marian, A. Kovner and H. Weigert, Phys. Rev.D59 014015 (1999); A. Kovner and J.G. Milhano, Phys. Rev. 
D61 014012 (2000) . A. Kovner, J.G. Milhano and H. Weigert, Phys.Rev. D62 114005 (2000); H. Weigert, Nucl.Phys. A 
703 (2002) 823. 

[12] E.Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A 692 (2001) 583; Phys. Lett. B 510 (2001) 133; E. Ferreiro, E. 

Iancu, A. Leonidov, L. McLerran; Nucl. Phys. A703 (2002) 489. 
[13] A. Mueller, Phys.Lett. B523 (2001) 243-248; 

[14] A. Mu eller and A.I. Shoshi; Nucl.Phys.B692: 175-208,2004; E. Iancu, A.H. Mueller and S. Munier; e-print archive 
hep-ph/0410018 

[15] E. Iancu and D. N. Triantafyllopoulos,e-print arxive hep-ph/0411405 
[16] A. Mueller, A. Shoshi and S. Wong; e-print arxive |hep -ph/0501088 
[17] E. Levin and M. Lublinsky; e-print archive hep-ph/0501173 

[18] J. Jalilian-Marian, A. Kovner, L. McLerran, and H. Weigert, Phys. Rev. D55 (1997) 5414-5428, 
[19] A. Kovner, M. Lublinsky and U. Wiedemann, in preparation; 
[20] A.Kovner and U.A. Wiedemann, Phys. Rev. D64:114002,2001; 

[21] A.M. Polyakov, Mod. Phys. Lett.A 3 (1988), 325; A.Yu. Alekseev and S.L.Shatashvili, Mod. Phys. Lett.A 3 (1988), 1551, 
D. Diakonov and V. Petrov; J.Exp.Theor.Phys.92:905-920,2001 e-Print Archive: hep-th/0008035 

[22] For the more detailed discussion of the relevant group theory see for example I. Kogan, A. Kovner and B. Tekin JHEP 
0105:062,2001; e-Print Archive: hep-th/0104047 

[23] To match the dimensions of t and x~ one should scale x~ by the longitudinal extent of the target. 

[24] See for example, Geometric phases in physics, A. Shapere and F. Wilczek editors, World scientific, Singapore, 1989 



